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Abstract 

Inspired by our previous finding that super symmetric Yang-Mills- Chern— 
Simons (SYM-CS) theory dimensionally reduced to 1+1 dimensions possesses 
approximate Bogomol'nyi-Prasad-Sommerfield (BPS) states, we study the anal- 
ogous phenomenon in the three-dimensional theory. Approximate BPS states 
in two dimensions have masses which are nearly independent of the Yang-Mills 
coupling and proportional to their average number of partons. These states are 
a reflection of the exactly massless BPS states of the underlying pure SYM the- 
ory. In three dimensions we find that this mechanism leads to anomalously light 
bound states. While the mass scale is still proportional to the average number of 
partons times the square of the CS coupling, the average number of partons in 
these bound states changes with the Yang-Mills coupling. Therefore, the masses 
of these states are not independent of the coupling. Our numerical calculations 
are done using supersymmetric discrete light-cone quantization (SDLCQ). 



1 Introduction 



It is extremely interesting to look for mechanisms in quantum field theory that give 
rise to anomalously light bound states, since these are the states we first see as we 
investigate new phenomena. Of course, spontaneous chiral symmetry breaking is the 
most profound example of these mechanisms. The underlying chiral symmetry of QCD 
gives rise to exactly massless bound states, the Goldstone bosons, which acquire masses 
from small mass terms for the constituent quarks which in turn break the symmetry. 

What we are suggesting here is that there may be a similar mechanism at work in 
supersymmetric theories. M = 1 pure supersymmetric Yang-Mills (SYM) theory in 
1+1 and 2+1 dimensions has massless BPS states ]1|, |2|. Recently || £|] we created 
an effective mass for the constituents of a SYM theory in 1+1 dimensions by adding 
a dimensionally reduced Chern-Simons (CS) term |5| to the theory. Although such a 
term does not break the supersymmetry, the theory with a CS term does not have BPS 
states. We found, however, that this theory has approximate BPS states, i.e. states 
whose masses are nearly independent of the Yang-Mills (YM) coupling. Furthermore, 
at large YM coupling the masses of these states are equal to the square of the CS 
coupling times the average number of particles. 

In this paper we will investigate the analog of these states in the (2+l)-dimensional 
theory. We will show that in 2+1 dimensions these states persist, although with some 
significant changes. We will see that in 2+1 dimensions the masses of these states are 
not independent of the YM coupling as they are in 1+1 dimensions. This is primarily 
because the average number of particles in these states grows with the coupling. How- 
ever, the masses of these states do not grow as fast as the masses of the other bound 
states, and, therefore, at strong coupling these states are anomalously light. 

The numerical method we will use is SDLCQ (Supersymmetric Discrete Light-Cone 
Quantization). This is a numerical method that can be used to solve any theory with 
enough supersymmetry to be finite. The central point of this method is that using 
DLCQ || [jj we can construct a finite dimensional representation of the superalge- 
bra ||. From this representation of the superalgebra, we construct a finite-dimensional 
Hamiltonian which we diagonalize numerically. We repeat the process for larger and 
larger representations and extrapolate the solutions to the continuum. 

In constructing the discrete approximation we drop the longitudinal zero-momentum 
mode. For some discussion of dynamical and constrained zero modes, see the review |7[] 
and previous work |J . Inclusion of these modes would be ideal, but the techniques re- 
quired to include them in a numerical calculation have proved to be difficult to develop, 
particularly because of nonlinearities. For DLCQ calculations that can be compared 
with exact solutions, the exclusion of zero modes does not affect the massive spec- 
trum 0. In scalar theories it has been known for some time that constrained zero 
modes can give rise to dynamical symmetry breaking M, and work continues on the 
role of zero modes and near zero modes in these theories ||10|| . Dropping the zero modes 
in CS theory does result in the loss of many of the interesting aspects of CS theory, 
most notably the quantization of the CS coupling. 

One particularly interesting property that will be preserved is the fact that the CS 
term simulates a mass for the theory. This property is interesting because reduced 
M = 1 SYM theories are very stringy theories. The low-mass states are dominated by 
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Fock states with many constituents, and as the size of the superalgebraic representation 
is increased, states with lower masses and more constituents appear [p], |], H [11], 0' 



b|, |TJ, [U| [IT], |T^, [LJ! pD| . The connection between string theory and supersymmetric 
gauge theory leads one to expect this type of behavior. However, these gauge theories 
are not very QCD-like, and ultimately one might like to make a connection with the 
low-mass spectrum observed in nature. The effective mass supplied by the CS term 
makes the theory more like QCD. 

The present Hamiltonian formalism provides us with the wave functions from which 
we construct the structure functions of the mass eigenstates. We will see that the 
structure functions of the approximate BPS states have unique shapes which might 
help to single them out experimentally. 

The paper is structured as follows. In Sec. |2| we provide a summary of SYM-CS 
theory; much of this is given in || but is repeated here for completeness. Our numerical 
techniques are described in Sec. An analysis of what should be expected for SYM- 
CS eigenstates is given in Sec. |], and our results are presented and discussed in Sec. ||. 
Section ^| contains a summary of our conclusions and some discussion of interesting 
further investigations. 

2 Supersymmetric Chern— Simons theory 

The Lagrangian of (2+l)-dimensional M = 1 SYM-CS theory is 

£ = Tr(-~£ Y M + i£ F + ~£cs), (1) 

where n is the CS coupling and 

£ym = F ftv F llv , (2) 
£ F = Vj^D^, (3) 

£cs = (a&Ax + jgApA v A^ + 2M . (4) 

The two components of the spinor \1/ = 2 _1//4 (^) are in the adjoint representation of 
U(N C ) or SU(N C ). We will work in the large- N c limit. The field strength and the 
covariant derivative are 

F ia , = d ll A v -d v A ll + ig[A ll ,A v ], D^ = d^ + ig[A^ ]. (5) 

The supersymmetric variations of the fields leads to the supercurrent in the usual 
manner via 

6C = ed li QM. (6) 

Light-cone coordinates in 2+1 dimensions are (x + ,x~,x ± ) where x + = x_ is the light- 
cone time and x x = —x±. The totally anti-symmetric tensor is defined by e + ~ 2 = — 1. 
The variations of the three parts of the Lagrangian in Eq. (|I]) determine the ('chiral') 
components Q ± of the supercharge via Eq. (H) to be 



/ d 2 xQ^ =(^ Q -)=~\j d 2 xe^ lxl + VF aP . 



(7) 
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Explicitly they are 

Q~ = -i2 i/A j d 2 x^(d + A- -d~A + + ig[A + ,A~ 

Q+ = -i2 5 ? 4 J d 2 xij(d + A ± -d 2 A + + ig[A + ,A ± 

One can convince oneself by calculating the energy-momentum tensor T^ v that the 
supercharge fulfills the supersymmetry algebra 

{g ± ,g ± } = 2 v / 2P ± , {Q + ,Q-} = -AP ± . (9) 

In order to express the supercharge in terms of the physical degrees of freedom, we 
have to use equations of motion, some of which are constraint equations. The equations 
of motion for the gauge fields are 

- D v F va = ^e auX F uX + 2g^ a V . (10) 

For a = + this is a constraint for A~ . In light-cone gauge, A + = 0, the constraint 
reduces to 

D_A- = ^-[(«- D 2 )D^ + 2(7*7+*]. (H) 

The equation of motion for the fermion is 7 M _D A1 \I / = — m*. Expressing everything in 
terms of ip and x leads to the equations of motion 

V2D + iP = (D 2 + k) X , (12) 
V2D_x = (D 2 -K)ip, (13) 

the second of which is a constraint equation. The constraint equations, (|TT| ) and (|TB|), 
are used to eliminate A~ and x- 

We compactify in the x 1 - direction with compactification length L, meaning that 
transverse momentum modes are summed over a discrete set of values 27m- 1 / L. In order 
to have a finite matrix representation for the eigenvalue problem, we must truncate 
these sums at some fixed integers ±T. The value of T defines a physical transverse 
cutoff = 2ttT/L; however, given this definition, T can be viewed as a measure of 
transverse resolution at fixed Aj_. 

At large N c the expansions of the field operators <fi = A± and ip, in terms of creation 
and annihilation operators for the Fock basis, are 

<j) ij (0,X~,X±) = 

y/2-KL , ^Jo y/2k+ L 

v n = — oo v 

if)ij(0,x~,x±) = 



Si* 



; 27rn^ 



^•(fc + , n ^) e -i^-+^-x + 6 t.( jfe + ) n ^) e ifc+ - ■ l 



From the field (anti-)commutators one finds 

a ij (p + ,n_ L ),a^(g + ,m_ L )] = {%(p + , nj_), b\ k {q + , m x )} = 5(p + - g + )5 n±im± ^5 jfc . (14) 
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The truncated supercharge Q can be written as 

Q~ = 9Qsym(T) + Q±(T) + mQ5 s (T) 

where 



Q1(T) 



2 3/4 



7T1 



\nM<T 



E j n dk^j= aljikjU^bijik,^) - b\ j (k,n ± )a ij (k,n ± ) 



and 



Qsym(^) 



| fc 2 - fci 
k 3 Vhk 2 

| fcl + ^3 

k 2 Vhh 

| ^2 + ^3 

kWhh 



|n x |<T 
£-5/4 



QcsCO = 2~ 1/4 i E /„ *^[at i (A:,n ± )6i i (A ; ,n ± ) + &t.(A ; ,n ± )a ii (A ; ,n ± ) 



(15) 

(16) 
(17) 



, 



/■CXD [ / 1 I 1 

E / dkxdk 2 dk 3 5(k 1 + k 2 - h)5 n ± +n ± n ± <2 — + — - — 

nM<T 1 2 3 

w r) & lj(^2, n^)b i:j (k 3 , rig ) + foJi( fc 3, )6 ifc (fci, n±)b kj (k 2 , n^)) 
4fc( fc i> n i) a lj( k 2i n 2 )hj{.h, jig ) - &2,.(fc3, ng )a ifc (A; 1 , n^)a kj {k 2 , n 2 ) 
a lk( k 3, n^)a kj (k u n^)b i:j (k 2 , n^) - a\ k (k u n^&j^Afe, n^a^k^ nj) 
b lk( k ii n i) a lj( k 2i n^)a i:j (k 3 , n^) - a^fa, nj)b ik (kx, n^)a hj {k 2 , n 2 ) 



(18) 



The symmetric truncation with respect to positive and negative modes aids in retaining 
a reflection parity symmetry in the states. The light-cone energy is (k\ + m 2 )/k + , and 
k± behaves like a mass. Here we see that k appears in a very similar way to k± and 
therefore behaves in many ways like a mass. 

When comparing the _L and CS contributions to the supercharge, we see that we 
have a relative i between them. Thus the usual eigenvalue problem 

2P+P» = V2P + (Q-) 2 \v) = V2P + (gQ SYM + Qj + inQ^fW) = M 2 n \ V ) (19) 

has to be solved by using fully complex methods. 

We retain^] the S'-symmetry, which is associated with the orientation of the large- N c 
string of partons in a state [21|. In a (l+l)-dimensional model this orientation parity 
is usually referred as a Z 2 symmetry, and we will follow that here. It gives a sign when 
the color indices are permuted 



Z 2 : aij(k,n J 



bij(k,n J 



-bji(k,n J 



(20) 



We will use this symmetry to reduce the Hamiltonian matrix size and hence the numer- 
ical effort. All of our states will be labeled by the Z 2 sector in which they appear. We 
will not attempt to label the states by their normal parity; in light-cone coordinates 
this is only an approximate symmetry. Such a labeling could be done in an approximate 
way, as was shown by Hornbostel fl2"2"| , and might be useful for comparison purposes if 
at some point there are results from lattice simulations of the present theory. 

1 We note that the CS term breaks transverse parity. 
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3 Numerical Methods 



We convert the mass eigenvalue problem 2P + P~\M) = M 2 \M) to a matrix eigenvalue 
problem by introducing a discrete P~ in a basis where P + and P 1 - = are diag- 
onal. As discussed in the Introduction, this is done in SDLCQ by first discretizing 
the supercharge Q~ and then constructing P~ from the square of the supercharge: 
P~ = (<5 _ ) 2 /v / 2- We have already introduced a finite discretization in the transverse 
direction, characterized by the compactification scale L and cutoff or resolution T. To 
complete the discretization of the supercharge, we introduce discrete longitudinal mo- 
menta k + as fractions nP + /K of the total longitudinal momentum P + . Here n < K 
and K are positive integers. The number of partons is also bounded by K. The integer 
K determines the resolution of the longitudinal discretization and is known in DLCQ 
as the harmonic resolution ||. The remaining integrals in Q~ are approximated by a 
trapezoidal form. The continuum limit in the longitudinal direction is then recovered 
by taking the limit K —>■ oo. In constructing the discrete approximation we drop the 
longitudinal zero-momentum mode. 

Our earliest SDLCQ calculations [lf§ were done using a code written for Mathe- 



matica and performed on a PC. This code was rewritten in C++ for the work pre- 
sented in [|14] and has been substantially revised to reach higher resolutions Q. The 



Mathematica code continues to be useful for checking any changes in the C++ code. 

To obtain the spectrum of the CS theory we solve the complex eigenvalue problem, 
Eq. flT9p . For the numerical evaluation we can exploit the structure of the supercharge 

where A and B are real matrices. The Hamiltonian has thus a simple decomposition 
into real and imaginary parts in the bosonic sector 



^boson = AA T + BB T + i (BA T - AB T ) , (22) 
and the fermionic sector 

^fcrmion = A T A + B T B + % (a t B - B T a) . (23) 



We extract several of the lowest eigenstates by applying the Lanczos algorithm [|23 
discussed in n. 



as 



4 Setting the stage 

To analyze the full SYM-CS theory, it is convenient to first study three subsets of the 
(2+l)-dimensional theory: the (l+l)-dimensional pure SYM theory (T — 0, k — 0), 
the pure SYM theory in three dimensions (T > 0, k = 0), and the (l + l)-dimensional 
SYM-CS theory (T = 0, k ^ 0). We have previously studied each of these theories 
separately 0, gjj [J |. 

Let us start with the dimensionally reduced pure SYM theory ||T9| , p0| , for which 



Q~ = 9Qsym(T = 0) . (24) 
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There are two unique properties that characterize this theory. As we increase the 
resolution, we find that there are new lower-mass states that appear, and this sequence 
of states appears to accumulate at M 2 = 0. In addition, there are massless BPS states. 
The dominant component of the wave function of the BPS states can be arranged 
to have 2, 3, . . . , K particles, i.e. up to the maximum number of particles allowed at 
resolution K. Therefore at resolution K there are K — 1 bosonic BPS states and K — 1 
fermionic BPS states. 

Next we focus on the case T > 0, k = 0, which is the pure SYM theory in 2 + 1 
dimensions 0. The supercharge Q~ becomes 

Q- = gQ SYM (T)+Ql(T). (25) 

Again there are a number of unique properties that characterize the theory. At small 
g and low energy we recover the (1+1) -dimensional theory. At higher energy we find 
a series of Kaluza-Klein states, which we discuss in detail elsewhere |]24||. At large 



coupling one might expect that the transverse Q±/ g term would freeze out, and one 
would see states that are a reflection of the (l+l)-dimensional theory. While we see 
these states, they are only a small part of the spectrum. Q$ym by itself wants to induce 
a large number of particles, as we discussed above. The contributions from gQ^yu an d 
are therefore minimized by a larger number of particles, each with a small transverse 
momentum. We therefore find that the average number of particles in the bound states 
grows with the coupling. The massless BPS states of the (l+l)-dimensional SYM 
theory persist in 2+1 dimensions, but now the number of particles in all of these states 



increases rapidly with the coupling. In Fig. 3 of Ref. (16j we displayed the average 
number of particles in the massless states as a function of the scaled coupling. At zero 
coupling we have the (l+l)-dimensional theory, and the states arrange themselves to 
have 2 through K particles. 

If we set T = we get SYM-CS theory reduced to two dimensions || . Here we find 
that the most important role of the CS term is to provide a mass for the constituents. 
This freezes out the long, lower-mass states that characterize (l+l)-dimensional SYM 
theory. Interestingly, however, the massless BPS states become massive approximate 
BPS states and have masses that are nearly independent of the YM coupling [|J, as 
seen in Fig. |l](a) where we show the two lowest Z2-even approximate BPS states have 
squared masses of 4k 2 and 16k 2 at infinite YM coupling. 

5 Numerical Results 



5.1 DLCQ analysis 

We start the discussion of our new results with a standard DLCQ analysis of the bound 
states at0 g = 0.5. At each value of the resolution K and in each Z2 sector we look 
at the ten lowest-energy bound states at each transverse resolution T. In order to 
extrapolate to infinite transverse momentum cutoff, we fit the curves of mass squared 
as a function of the inverse transverse cutoff with a linear function in 1 /T, for T > 3, as 
shown in Fig. 0(a). The intercepts of these curves are the transverse continuum results, 
which we plot as a function of 1/K in Fig. 0(b). Again we make a fit to the lowest 

2 Values of g are given in units of y/4ir 3 /N c L. 
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0.6 0.8 



(a) (b) 

Figure 1: The Z 2 -even spectrum of SYM-CS theory as a function of the Yang-Mills 
coupling g in (a) 1+1 dimensions and (b) 2+1 dimensions. For (a), the longitudinal 
resolution is K — 9, M 2 is measured in units of k 2 , and g is in units of k^Jti/N c . 
For (b), the longitudinal resolution is K — 6, the transverse resolution is T = 1, the 
Chern-Simons coupling is k = 2tt/L, M 2 is in units of 4rr 2 /L 2 , and g is in units of 
An 3 /N C L. 




0.2 0.4 0.6 0.8 1 . 0.1 0.2 0.3 

1/T 1/K 

(a) (b) 

Figure 2: The i^-even spectrum of (2+l)-dimensional SYM-CS theory as a function 
of (a) transverse resolution T and (b) longitudinal resolution K. The values in (b) 
have been extrapolated to T = oo. The coupling strengths are g = 0.5, in units of 
^4tt 3 /N c L, and k = 2n/L. 



mass state for K > 4, this time linear in 1/K, and the intercept of this curve is the 
true continuum mass. We have not attempted to fit the 1/K behavior of the higher 
states, because we are mainly interested in the lowest state, which is the analog of the 
lowest approximate BPS state in two dimensions. 

For the case g = 0.1 the average number of particles is well behaved, in the sense 
that it stays small. This allows for an additional approximation. Namely, we can 
restrict the number of partons per Fock state. This is a good approximation to the 
lowest ten states. Their average parton number is either two or three for even or odd 
Z2, respectively. This is clear since the states at small coupling are very close to the 
free states, which in turn are the discrete manifestations of the two and three particle 
continua starting at M 2 = 4k 2 and M 2 = 9k 2 , where k 2 is the effective parton mass. 
We restrict the number of partons in a state to four. This additional approximation 
allows us to carry the calculation to nine units of transverse resolution, T = 9, and nine 
units of longitudinal resolution, K = 9. The details of this analysis will be presented 
elsewhere |[24j| . At g = 0.5 we are unable to make this approximation and are only able 



to go to maximum resolutions of K = 6 and T = 4. The reason that this approximation 
breaks down is that, as we increase T and K, the average number of partons in a state 
increases substantially. When this number approaches K, we miss a significant part of 
the wave function in the calculation. 

In the Z2 = — 1 bosonic sector the lowest state is a three-particle state at weak 
coupling. All the analysis of this sector goes through similarly to the Z2 = +1 sector 
and will also be presented elsewhere [El . 
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5.2 Dependence of the spectrum on the Yang-Mills coupling 

In Fig. [lj(b) we show the coupling dependence of SYM-CS theory in 2+1 dimensions 
in the Z 2 = +1 sector. We will restrict ourselves to this sector in the discussion that 
follows. For numerical convenience the results presented here are at specific and low 
values of the longitudinal and transverse cutoffs K and T but at many values of the 
coupling. The point we want to make is that the phenomenon of approximate BPS 
states persists from two dimensions: some states appear to be nearly independent of 
the coupling at fixed T and K. 

We can distinguish two regions in the low-lying spectrum. In the region of very 
small g we have a spectrum of almost free particles, and in the regime of larger g an 
approximate BPS state is clearly visible while the masses squared of the rest of the 
states grow like g 2 . In the intermediate region we see a large number of level crossings 
as the states rearrange themselves. This transition occurs around g = 0.3. We will 
refer to couplings above this cross-over point as being strong, because the spectrum 
has completely rearranged itself here. 

We have investigated the spectrum at g — 0.1 and g = 0.5 as a function of K and 
T as typical representatives of weak and strong coupling, respectively. It is obvious 
that some of the states seen in Fig. |l|(b) move differently as a function of K and T. At 
weak coupling we have the manifestation of a multi-particle continuum within a discrete 
approach, and the masses of the lowest states are very close together. As the coupling 
grows, a large gap develops between the lowest-mass state, which is an approximate 
BPS state, and the other states in the spectrum, as seen in Fig. |](b). The anomalous 
behavior of the next to lowest state, which is the only state to decrease in mass, is a 
small tranverse cutoff effect. When T is increased the gap is restored, as can be seen in 
Fig. |K a): the mass of the next to lowest state at T = 1 increases rapidly as T grows, 
leaving a mass gap in the transverse continuum limit. The computing time required to 
produce a graph of the spectrum as a function of the coupling at large T is prohibitive, 
but we display the lowest masses in Table [I] at g = 0.1 and g = 0.5. We see that at 
g = 0.1 the lowest 10 states are more or less uniformly spaced from M 2 = 4 to M 2 = 10 
in units of 4tt 2 /L 2 . This is what we would expect since at small coupling these states 
are closely related to the two-particle continuum states. At strong coupling, g = 0.5, 
we see that a large gap has developed between the state at M 2 = 8.93 and M 2 = 13.14. 
These results are obtained at resolution K = 6, but from Fig. 0(b) it appears that this 
gap is growing with the longitudinal resolution, and, therefore, in the continuum we 
expect it will be even larger. 

We have studied the average number of particles of all the states in all of the sectors 
and a detailed discussion will be presented elsewhere [02] . Here it is interesting to note 
that for the state at M 2 = 8.93 we find an average parton number squared (n 2 ) = 8.23. 
This is close to the expected value based on the duality discussed in SYM-CS in 1+1 
dimensions 0] , where we found that the mass at infinite coupling is the average parton 
number squared at CS coupling unity. By contrast, the average parton number squared 
(n 2 ) = 9.9 of the second lightest state is far away from its mass squared M 2 = 13.14. 
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M' 2 


0.1 


4.36 


4.85 


4.93 


7.41 


7.47 


8.35 


8.36 


9.34 


9.43 


9.48 


0.5 


8.93 


13.14 


13.88 


14.92 


15.10 


15.62 


15.82 


17.38 


18.35 


19.76 



Table 1: Mass squared in units of 4tt 2 /L 2 of the lowest-mass states at g = 0.1 and 
g = 0.5, in units of ^J4tt 3 /N c L, showing the development of a gap above the lowest 
state. Results in both cases are for resolutions K — 6 and T = 4 and for Chern-Simons 
coupling k = 2n/L. 

5.3 Structure Functions 

In this section we will present a discussion of the structure functions of the approximate 
BPS states. A complete discussion of the structure functions of all the states will be 
presented elsewhere ||24| . The structure functions are functions of both the longitudinal 
and transverse momentum and are defined to be the probability that a particle of type 
A will have a longitudinal momentum fraction x = k + /P + and transverse momentum 
k- 1 . It is given as follows in terms of the light-cone wave function ip: 

g A {x,k L ) = dx! ■ ■ ■ dx q J dk^ ■ ■ ■ dk^d (^1 x * ~ ^ S H 

x K*i ~ ~ |^(x ls ki; . . . x q , k£) | 2 . (26) 

i=i 

We present the discrete version of the structure function as a function of the number 
of units of longitudinal momentum and the number of units of transverse momentum. 
No attempt has been made to extrapolate these curves. 

In Fig. |3] we see the structure function of the approximate BPS state of the SYM- 
CS theory for g = 0.1 and g = 0.5. At weak coupling it has one localized peak in 
both variables. In fact this is true of all of the low-mass states of this theory at weak 
coupling, and the approximate BPS state does not appear in any way special. As we 
move to stronger coupling, however, all structure functions of the normal bound states 
remain peaked in nil, while the peak moves down to lower values of n\\ as the number of 
particles increases. In contrast we see that the BPS state is nearly flat in Tin, although 
it has a small peak at n\\ = 1. This behavior is unique to this state. 

At this point we do not have an understanding of why the approximate BPS state 
has this behavior. It is, however, clearly important that it does, if for no other reason 
than to give one a means to identify such a state experimentally. 

6 Outlook 

In our earlier work on dimensionally reduced SYM-CS theory []3], Q , we found that there 
were states whose properties are largely dictated by the properties of the massless BPS 
states of the underlying SYM theory. We saw that the masses of these states were 
nearly constant as functions of the YM coupling, and we therefore called these states 
approximate BPS states. The mass scale of these states is set by the CS coupling and 
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(a) (b) 

Figure 3: Structure function of the lowest approximate BPS state at (a) g = 0.1 and 
(b) g = 0.5, in units of ^JAtt 3 /N c L. The numerical resolutions are K = 6 and T = 4. 
The value of the Chern-Simons coupling is k = 2n /L. 

the average number of partons. The binding energy is approximately zero, because of 
the underlying properties of the BPS state. 

In this paper we have looked for the effect of these approximate BPS states in the 
full (2+l)-dimensional SYM-CS theory. We have looked at the low mass spectrum 
and found that at strong coupling there are states with an anomalously low mass. The 
mass scale for these states is set by the CS coupling and the average number of partons, 
similar to the (l+l)-dimensional case. These states are a reflection of the exact BPS 
states that exist in (2+l)-dimensional SYM theory. In addition, we examined the 
structure functions for these states. We find that they are different from the structure 
functions of the other states; they are nearly flat in longitudinal momentum. 

It appears that we have found a new mechanism for generating low-mass bound 
states. There are two issues, however, that need to be addressed before speculating 
on the possibility that this mechanism might be realized in nature. First, we have not 
included any matter fields in this theory. There has been some work done with matter 
fields in SDLCQ [PH| , and one of our next projects will be to include them. Second, we 
have not included supersymmetry breaking. Instead, we added masses for the adjoint 
fields by including a CS term. The calculational technique we use, SDLCQ, relies 
heavily on the fact that we have an exactly supersymmetric theory, and it is essential, 
therefore, that we introduce masses without breaking the supersymmetry. Inclusion of 
supersymmetry breaking remains an outstanding problem in this numerical approach. 
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